A nonsmooth Newton multigrid method for a hybrid, shallow model of marine ice sheets
نویسندگان
چکیده
The time evolution of ice sheets and ice shelves is modelled by combining a shallow lubrication approximation for shear deformation with the shallow shelf approximation for basal sliding, along with the mass conservation principle. At each time step two p-Laplace problems and one transport problem are solved. Both p-Laplace problems are formulated as minimisation problems. They are approximated by a finite element truncated nonsmooth Newton multigrid method. As an illustration, we compute the steady state shape of an idealized ice sheet/shelf system.
منابع مشابه
An adaptive Newton multigrid method for a model of marine ice sheets
Article history: Received 5 December 2012 Received in revised form 26 May 2013 Accepted 25 June 2013 Available online 4 July 2013
متن کاملSeminar of Numerical Analysis
For more than one decade, Antarctic ice sheets have been retreating dramatically, and are expected to shrink even more quickly in the future. Recent studies have identifed the migration of the grounding lines separating the grounded part of ice sheets and the floating part of ice shelves as being a key process controlling marine ice sheet stability [4]. It is also wellknown that the grounding l...
متن کاملAchieving Textbook Multigrid Efficiency for Hydrostatic Ice Sheet Flow
The hydrostatic equations for ice sheet flow offer improved fidelity compared with the shallow ice approximation and shallow stream approximation popular in today’s ice sheet models. Nevertheless, they present a serious bottleneck because they require the solution of a 3D nonlinear system, as opposed to the 2D system present in the shallow stream approximation. This 3D system is posed on high-a...
متن کاملSolution of nonlinear Stokes equations discretized by high-order finite elements on nonconforming and anisotropic meshes, with application to ice sheet dynamics
Motivated by the need for efficient and accurate simulation of the dynamics of the polar ice sheets, we design high-order finite element discretizations and scalable solvers for the solution of nonlinear incompressible Stokes equations. In particular, we focus on power-law, shear thinning rheologies commonly used in modeling ice dynamics and other geophysical flows. We use nonconforming hexahed...
متن کاملTruncated Nonsmooth Newton Multigrid Methods for Simplex-Constrained Minimization Problems
We present a multigrid method for the minimization of strongly convex functionals defined on a finite product of simplices. Such problems result, for example, from the discretization of multi-component phase-field problems. Our algorithm is globally convergent, requires no regularization parameters, and achieves multigrid convergence rates. We present numerical results for the vector-valued All...
متن کامل